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p ,' Abstract 

Can dark matter be stabilized by charge conservation, just as the electron is in the standard 

I model? We examine the possibility that dark matter is hidden, that is, neutral under all standard 

' model gauge interactions, but charged under an exact U(l) gauge symmetry of the hidden sector. 

' Such candidates are predicted in WIMPless models, supersymmetric models in which hidden dark 

. matter has the desired thermal relic density for a wide range of masses. Hidden charged dark 

jy-^ ! matter has many novel properties not shared by neutral dark matter: (1) bound state formation 

' and Sommerfeld-enhanced annihilation after chemical freeze out may reduce its relic density, (2) 

ON ' 

o 
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similar effects greatly enhance dark matter annihilation in protohalos at redshifts of 2 ~ 30, (3) 
Compton scattering off hidden photons delays kinetic decoupling, suppressing small scale structure, 
and (4) Rutherford scattering makes such dark matter self-interacting and collisional, potentially 
impacting properties of the Bullet Cluster and the observed morphology of galactic halos. We 
I analyze all of these effects in a WIMPless model in which the hidden sector is a simplified version 

of the minimal supersymmetric standard model and the dark matter is a hidden sector stau. We 
find that charged hidden dark matter is viable and consistent with the correct relic density for 
reasonable model parameters and dark matter masses in the range 1 GeV ^ mx ^ 10 TeV. At the 
same time, in the preferred range of parameters, this model predicts cores in the dark matter halos 
of small galaxies and other halo properties that may be within the reach of future observations. 
These models therefore provide a viable and well-motivated framework for collisional dark matter 
with Sommerfeld enhancement, with novel implications for astrophysics and dark matter searches. 

PACS numbers: 95.35.-hd, 04.65.+e, 12.60.Jv 
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I. INTRODUCTION 



Dark matter makes up most of the matter in the Universe, but its identity is unknown. 
So far, dark matter has been observed only through its gravitational interactions. A logical 
possibility, then, is that dark matter is hidden, that is, neutral under all standard model 
(SM) gaug e interactions. This possibility has been explored for many years H, 0, M, @, 
0, Sis, IG, m and brings with it a great deal of model building freedom. Unfortunately, this 
freedom is generically accompanied not only by a loss of predictivity, but also by the loss of 
appealing features common to other dark matter candidates, such as the strong connections 
to central problems in particle physics, such as the gauge hierarchy and flavor problems, and 
the naturally correct thermal relic density of weakly-interacting massive particles (WIMPs). 

In this work, we consider the possibility that hidden dark matter is charged under an 
exact U(l) gauge symmetry in the hidden sector. We assume the hidden sector is truly 
hidden, with no "connector" particles charged under both SM and hidden gauge groups. 
Although many of our results are relevant for any charged hidden dark matter, including 
cases with connectors or massive hidden photons, we are particularly motivated by the 

ITI . which ameliorates some of the 
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problems of hidden dark matter noted above. (For other recent discussions of hidden dark 
matter charged under exact gauge symmetries, see Refs. 0, H, iH.) In this framework, 
dark matter is hidden, but additional structure implies that it nevertheless naturally has the 
correct thermal relic density. In the examples studied in Refs. BQQEHI |l7j , the dark 
matter is in a hidden sector of supersymmetric models with gauge-mediated supersymmetry 
breaking (GMSB). Assuming that supersymmetry is broken in a single sector, the hidden 
sector mass scale mx and gauge couplings gx are related to the observable sector weak scale 
mw and gauge coupling by 

mx mw 
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9w 



This implies that the thermal relic density of WIMPless dark matter, Vt 
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/ g^ ~ f^wiMP, and so is naturally of the desired magnitude, even though mx 



need not be at the weak scale. The WIMPless framework also has other virtues, including 
that fact that it is naturally consistent with an elegant solution to the new physics flavor 
problem, since the gravitino mass, and with it all dangerous flavor- violating gravity-mediated 
effects, satisfies ttlq -C mw- This contrasts with the case of neutralino dark matter, where 
stability of the neutralino y requires ttiq > m^ ~ mw- Other GMSB dark matter candidates 



have been proposed [22|, |23|, |2J, |25!, |26 



However, it is notable that the relation mx oc g\ 
results directly from the desire to generate a flavor-blind superpartner spectrum. WIMPless 
dark matter therefore not only dissolves the tension between dark matter and the flavor 
problem in supersymmetry, but it exploits an elegant solution of the supersymmetric flavor 
problem to predict the correct thermal relic density. 

Assuming mQ < mx, however, how is the stability of the hidden sector candidate main- 
tained? Although there are many possibilities in the WIMPless framework [S], an elegant 
possibility is that the dark matter is charged under an exact gauge symmetry in the hid- 
den sector. The dark matter candidate may then be stable for the same reason that the 
electron is stable — its decay is prevented by charge conservation. We are therefore led 
rather straightforwardly to the possibility of charged dark matter in a hidden sector. In the 
following sections, we will consider, as a concrete example, a minimal supersymmetric stan- 
dard model (MSSM)-like hidden sector with C(l) Yukawa couplings. The massless hidden 
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particles — the hidden photon, gluon, and neutrinos — are then all neutral under hidden 
U(1)em- Any charged particle in the hidden sector, for example, the hidden W, chargino, 
stau, tau, or charged Higgs boson, would then be stable, provided it was the hghtest mx 
particle. A particularly plausible choice is the lighter stau f, whose mass receives smaller 
SUSY-breaking contributions than other hidden superparticles. 

This possibility of charged hidden dark matter requires a rethinking of many issues in 
cosmology, as its properties differ markedly from the more conventional possibility of neutral 
dark matter. In this work, we analyze these issues comprehensively to establish the viability 
of the scenario and determine what ranges of parameters are allowed. Because the dark 
matter is charged, it may form bound states and then annihilate in the early Universe. 
Dark matter annihilation is also enhanced by the Sommerfeld effect 28|. The formation 



of bound states and the Sommerfeld effect are possible also in the standard model; for 
example, bound state formation impacts e~^e~ annihilation to 511 keV photons at late times 
when the e"*" and e~ are non-relativistic. In the present case, both effects may reduce the 
thermal relic density, potentially destroying the nice relic density property of WIMPless 
dark matter. After chemical freeze out, the dark matter's charge also plays an important 
role in kinetic decoupling, since charged dark matter may be far more efficiently coupled to 
the thermal bath of hidden photons through Compton scattering than neutral dark matter. 
This impacts the matter power spectrum, requiring a modification of the standard WIMP 
analysis [2^, 30, 31|. This also changes the expected minimum mass of dark matter halos. 
It is important to check the minimum halo mass because the Sommerfeld enhancement 
charged dark matter annihilation is most effective in the least massive halos where relative 
velocities are the smallest, and this effect has been noted to provide a stringent constraint 
on the possibility of charged dark matter [3^. Finally, the existence of a long range force 
between dark matter particles implies dark matter is coUisional. It is therefore also critical 
to investigate whether this scenario is consistent with constraints on self-interacting dark 
matter from the Bullet Cluster [33], observed ellipticity of dark matter halos |3J] and bounds 
from considerations of galactic dynamics [l9| . 

In the end, we will find that all of these constraints are satisfied for reasonable model 
parameters, and are consistent with the correct relic density for dark matter masses in the 
range 1 GeV ^ mx ^ 10 TeV. These dark matter candidates are therefore viable, sharing 
many of the theoretical virtues of WIMPs, but predicting drastically different features for 
astrophysics. In its minimal form, this model also predicts no direct detection, indirect 
detection, and collider signals. If desired, however, with the addition of connector sectors, 
the charged hidden dark matter model we consider here may also have implications for 
direct searches, indirect searches, and particle colliders. Some of these implications have 
been considered previously 
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and others are possible [35|, l36|, l37|, l38l . 



ISL- 

50(1 . We will not expand on these possibilities here, 
other than to note that our results show that an exact hidden U(1)em is possible, leading to 
a simple framework in which annihilation rates are enhanced by the pure Sommerfeld effect. 

This paper is organized as follows. In Sec. [ITl we briefly review our hidden sector model 
and its parameters, and discuss its thermal relic abundance. In Sec. Illlt we study the 
kinetic decoupling of the hidden sector dark matter, and in Sec. IIVI we present the resulting 
matter power spectrum and derive the minimum mass of dark matter halos. Details of this 
derivation are given in Sec. \M The impact of bound states and enhanced annihilation on 
relic densities is evaluated in Sec. |Vl and the constraint from annihilation in protohalos is 



discussed in Sec. |Vll In Sec. IVIH we study the effect of the long range force between dark 
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matter particles on galaxy shapes and derive bounds on the parameter space. In Sec. IVIIII 
we verify that our model is consistent with bounds on self-interactions, such as those from 
merging clusters. We summarize our conclusions in Sec. IIXI We note that this work is 
necessarily rather wide-ranging; readers interested in only a subset of the topics may find 
Sec. [IX] a useful guide to the relevant conclusions and figures. 



II. HIDDEN SECTOR: PARTICLES, PARAMETERS, AND RELIC DENSITY 



In this work, following Ref. [16|, we consider a hidden sector that is MSSM-like, with 
gauge group SU(3) x S'U{2)l x U(l)y, but with only one generation of fermions and 0{1) 
Yukawa couplings. It is important that the model be chiral, so that hidden electroweak 
symmetry breaking sets the scale for the masses and the naturally correct thermal relic 
density is realized. Of course, the model should also be anomaly-free. The one-generation 
MSSM is a relatively simple model that satisfies these criteria, and has the advantage that 
one can apply much of the intuition built up from studies of the usual MSSM. Many other 
choices are possible, however, and would also be worth further study. Assuming electroweak 
symmetry breaking that leaves SU(3) x U(1)em preserved in the hidden sector, the particle 
masses in the hidden sector are 

Hidden weak scale : W'', Z^, , g'\ x^^, t^, b'', t^, h'^, f ^, 

~0: 7',/,^',G, (2) 

where we have used third-generation notation for the hidden (s)quarks and (s)leptons to 
remind us that they have 0{1) Yukawa couplings. From this spectrum, a natural possibility 
is that the lightest weak scale particle is charged, and therefore stabilized by U(1)em charge 
conservation in the hidden sector. As noted in Sec. [U there are many possible charged 
dark matter candidates in the hidden sector. To be concrete, we will focus here on the 
right-handed stau, which we denote from now on. 

Despite the complexity of the hidden sector, with a few mild assumptions, there are only 
a few parameters that are relevant for this study. These are 

mx, rrizh, ax, tan6'(V, ^rh = , (3) 

J-KH 

where mx = frifh is the dark matter particle's mass, m^h is the mass of the hidden Z boson, 
c^x = {e^Y /An is the fine structure constant of the hidden U(1)em, Ow the hidden sector's 
weak mixing angle, and ^rh is the ratio of hidden to observable sector temperatures at the 
time of reheating. The parameter ^rh need not be 1; for example, the infiaton may couple 
differently to the observable and hidden sectors, leading to a temperature asymmetry [1, 3] ■ 
For a one-generation MSSM to satisfy Big Bang nucleosynthesis constraints on massless 
degrees of freedom, one requires ^rh ^ 0.8 [16]. We will consider values in the range 
0.1 < ^RH < 0.8. In the hidden sector, the 7^^ and temperatures are always identical, and 
we denote the common temperature by T^. In the visible sector, the 7 and v temperatures 
diverge eventually, of course; we denote let T denote the 7 (cosmic microwave background 
(CMB)) temperature. We do not assume hidden sector gauge unification, although we 

will present some of our results for the grand unified value ianO'^ = -i/3/5 ~ 0.8 as an 



interesting example. We assume that the hidden sector's strong coupling gg is small enough 
that hadronization effects are negligible in the early Universe. 
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The hidden dark matter's rehc abundance is determined by the parameters of Eq. 



The dark matter annihilates through r^+r^- 
handed, its thermally-averaged annihilation cross section is |l6l 



If the r is dominantly right- 
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To avoid presenting results for the special case in which the annihilation cross section is 
resonantly enhanced by the pole, we fix = 1.5mx- The dependence of (cr^f) on 
Tzh is therefore negligible. In the numerical studies below, we will consider weak mixing 



angles in the range yS/S < ta,n9^r < 10. If the hidden sector has custodial symmetry, 
m-iyh = rrizh cos 9^, then tan 9'^ > 1.11 implies < nifh. Although would also be 
a perfectly good charged hidden dark matter candidate, we will instead assume custodial 
symmetry is broken, for example, by triplet Higgs bosons, so that f'^ remains the dark 
matter candidate for all values of tan^jy. 

With these assumptions, contours of constant Qxh"^ in the {mx,<yx) plane are given in 
Fig. [H (Constraints from observations are also shown in Fig. [T| these will be discussed in 
Sees. IVIII and IVIIII ) As evident from Fig. [1] and Eq. ([5]), there is a strong dependence on 
t3Jl9y^r^. for fixed ax, the neutrino cross section is enhanced by 1/ cos'^ 9^y for large tan^{y. 
This dependence is very important. In high energy processes, such as thermal freeze out, 
both EM and weak processes are effective. In low energy processes, such as those occurring 
at present, only EM interactions mediated by the massless 'y^ are important. If there were 
no hidden weak interactions, for fixed mx and ^rh, flxh"^ — 0.11 would fix ax, and so 
the relic density would fix the strength of dark matter self- interact ions now. In the more 
general case we are considering, however, one can keep fix constant for arbitrarily small 
ax, provided ta.n9^, is large enough. As noted above, we will consider tan6'(j^ as large as 
10; larger values are, of course, also possible. 



III. KINETIC DECOUPLING OF DARK MATTER 



In the standard cosmology, structure formation is hierarchical, that is, small objects 
form first and progressively the larger objects form via mergers and accretion. The mass 
of the smallest dark matter halo is determined by the dark matter particle's mass and the 
temperature at which it kinetically decouples from the other light particle species. For 
a typical WIMP with electroweak scale mass, the kinetic decoupling temperature ranges 
roughly from 10 MeV to 1 GeV, which leads to the formation of dark matter halos with 
masses as small as lO'^^M© to IO'^Mq |29|, [30|, [mI, [mI, [Hi, [H, Ei, [HI]. For other dark 
matter particles, such as weak-scale gravitinos or MeV dark matter, the minimum mass 
dark matter halo could be as large as the smallest dwarf galaxies 5^, 57, 5^. Some regions 
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FIG. 1: Allowed regions in {mx,(Xx) plane, where mx is the mass of the dark matter charged 
under the unbroken hidden sector U(1)em with fine-structure constant ax- Contours for fixed 
dark matter cosmological relic density consistent with WMAP results, ^xh'^ = 0.11, are shown 
for (tan0(^,^RH) = (a/3/5, 0.8), (v^3/5, 0.1), (10,0.1) (dashed), from top to bottom, as indicated. 
The shaded regions are disfavored by constraints from the Bullet Cluster observations on self- 
interactions (dark red) and the observed ellipticity of galactic dark matter halos (light yellow). 
The Bullet Cluster and ellipticity constraints are derived in Sees. IVIIII and IVIIl respectively. 



of the parameter space of these models are excluded because the predicted minimum mass 
halo is in conflict with observations. 

In this section, we analyze the kinetic decoupling of hidden charged dark matter. One 
notable difference between the WIMP and hidden charged dark matter is that the charged 
dark matter interacts not only through weak interactions, but also through EM interactions. 
For the case of dark matter, this implies that the dark matter remains in kinetic contact 
not only through the weak process f^z/'* ^ f'^u^, but also through the Compton scattering 



process r^'j^ 



f 7 . As we will see, at low temperatures, the thermally-averaged weak cross 
section is suppressed by T^^/m^, but this suppression is absent for Compton scattering, 
creating a large, qualitative difference between this case and the canonical WIMP scenario. 
Note also that, in principle, in the case of charged dark matter, bound state formation also 
impacts kinetic decoupUng. As we will see in Sec. |Vl however, very few staus actually bind, 
and so this effect is not sig nificant and may be neglected in our analysis. 
We follow Refs. 1 



55l | to determine the temperature of kinetic decoupling for the dark 



matter particle. In the hidden sector, the Boltzmann equation governing the evolution of 
the dark matter particle's phase space distribution is 



dm 

dt 



r(T'*)(T"mxAp^ + p-Vp-+3)/(p-) 



(6) 
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where 



is the momentum transfer rate. In Eq. ([7]), Qq is the number of degrees of freedom of the 
scattering particle {qq = 2 for both and i/^), the A''^_,_3 are constants defined in Ref. 55 
and the c„ are determined by parameterizing the colhsion amphtudes, evaluated at t 
and s = m\ + 2mxT^-, as 
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t=0,s=m\+2mxTt^ 
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(8) 



In the present case, with the help of the CalcHEP program we find that the 

squared amplitudes of the relevant processes are 



2 QAn^ax [(^3^ - s) 



2« 



syst + {m\ 



lQTi^a\ [{m\ - sf + St] 



cos^ 6'{y (m| 



and so 



Casting these results in the form of Eq. 



t=0,s=m\+2mxT'^ 
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\mx^ 



(9) 
(10) 

(11) 
(12) 



45 rrix 



, we find 

1247r5a^ mx T'^^ 
63 cos^ Oy^ m% 



mx, 



(13) 



The inverse of T{T^) is the time needed for the hidden sector photons and neutrinos 
to transfer momentum |p| ~ to the staus and thereby keep the dark matter in kinetic 
equilibrium. During this phase, the temperature of the dark matter tracks that of the 
hidden section photons or neutrinos, depending on the dominant scattering mechanism. 
The contribution from Compton scattering (the first term) does not suffer the {T^/mxY 
suppression typical of weak interactions. Therefore, for tan6'{!i/ ~ 1, Compton scattering 
dominates the contribution to T{T^). However, for tan6'{y = 10, the channel is enhanced 
dramatically and dominates for mx ^ 1 GeV. 

We define the kinetic decoupling temperature as the temperature where the momentum 
transfer rate and the Hubble expansion rate become equal, that is, T{T^^) = H{T^^). (For 
other possible definitions, see, for example, Refs. [5J, l55|.) Here we adopt the hidden sector 
point of view and write the Hubble parameter as 



/47r3 
I5" 



T 



h2 



pi 



(14) 



where gl°^{T^) = g^^ + g^, (t/T^^ is the total effective relativistic degrees of freedom, and T 
is the visible sector's photon temperature. We assume there is no thermal contact between 
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FIG. 2: Kinetic decoupling temperatures as functions of dark matter mass mx for (tan 6yy, ^rh) = 
(10,0.1), (^^3/5, 0.1), (1/3/5, 0.8), from top to bottom, as indicated. For each case we plot both 
the hidden sector photon temperature T^^ (dashed) and the corresponding visible sector photon 
(CMB) temperature Tkd (solid) at the time of kinetic decoupling. 

hidden and observable sectors, and so the visible sector enters the analysis only through the 
contribution of its effective degrees of freedom to the expansion rate. 

In Fig. [2] we show the hidden sector kinetic decoupling temperature T^^ and visible sec- 
tor kinetic decoupling temperature T^d as functions of mx for various combinations of 
(tangly, ^rh)- The hidden U(1)em fine-structure constant ax is determined by requiring 
the correct dark matter relic abundance Vtxh"^ = 0.11. To understand these results, consider 

first the case (tan6'(j(/, ^rh) = (-\/3/5, 0.8) and nix ~ 100 GeV. In this case, gi ~ (?2 in the 
hidden sector, the hidden and observable sectors have comparable temperatures, and the 
dark matter has a weak-scale mass. This case is thus similar to the standard WIMP case, 
except that the hidden dark matter is charged. For the WIMP case, one typically finds 
Tkd ~ 10 MeV — 1 GeV, as noted above, whereas in this charged hidden dark matter case, 
we find that Tkd ~ 0.1 MeV. We see that the presence of Compton scattering does in fact 
have a large impact, keeping the charged dark matter in kinetic equilibrium to much lower 
temperatures. 

This effect is moderated for other values of (tan 6^r, ^rh)- For (tan 6'{^, ^rh) = (^3/5, 0.1), 
the hidden sector is colder relative to the visible sector. As shown in Ref. [3], a colder hidden 
sector requires a smaller ax to meet the relic abundance constraint, which results in a smaller 
momentum transfer rate. Furthermore, for a colder hidden sector, the visible sector appears 
hotter. The resulting larger Hubble expansion rate of Eq. also causes the momentum 
transfer processes to become inefficient earlier. This effect is even more pronounced for 
(tan^l?!^, ^rh) = (10,0.1), since large taxiO^ implies even lower ax, increasing the visible 
sector Tkd further. 
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IV. POWER SPECTRUM OF DARK MATTER DENSITY FLUCTUATIONS 



The interactions and subsequent decoupling of dark matter particles damp the matter 
power spectrum. The matter power spectrum is the Fourier transform of the two-point 
function of the cosmological density fluctuations in matter. For Gaussian initial conditions, 
this encodes all the relevant information for linear perturbation theory. In this section, we 
will study the power spectrum of hidden charged dark matter to contrast its predictions with 
that of canonical WIMP models. We will focus on small scales where differences might be 
expected, especially in the predictions for the mass of the smallest dark matter structures 
that can form. We will find that the differences from the usual WIMP models are not 
observable with current data, although these may be observable in the future. However, 
the smallest halo mass prediction is important for the calculation of Sommerfeld-enhanced 
self-annihilations, which we will discuss in the following section. 

The damping due to the coupling of charged dark matter to hidden sector photons (and 
neutrinos) results from two distinct processes. The coupling of the dark matter to the 
hidden sector photons and neutrinos leads to damped oscillatory features [52] due to the 
interplay between the gravity of dark matter and the pressure of the photons. This is the 
dominant effect for the case of WIMPs with electroweak scale masses and also for our case 
with mx > 1 GeV. In addition, after decoupling, the free-streaming of the dark matter 
particles further suppresses the power spectrum. This effect dominates at lower masses and 
in particular in the region of parameter space where interactions with the hidden sector 
neutrinos dominate over those with hidden sector photons. 

The presence of the hidden sector and the coupling of the dark matter to 7^ and 
make the decoupling and matter power spectrum calculation different from the standard 
WIMP case. The two main differences are the following. In the standard WIMP case, the 
WIMP couples to standard model fermions that are part of a tightly coupled (coUisionally- 
coupled) fluid. This implies that multipole moments of the density fluctuations higher than 
the dipole are strongly suppressed. This result stems from the fact that in a coUisionally- 
coupled plasma, a quadrupole anisotropy pattern can only develop if there is significant 
diffusion. In the present case, the hidden sector photons and neutrinos are only coupled to 
the dark matter and are able to diffuse significantly. We therefore need to track their density 
fluctuations to higher multipoles. The other notable difference is that the expansion rate is 
predominantly set by the visible sector matter and so the ratio of the scattering rate in the 
hidden sector to expansion rate is different from the usual WIMP case. We provide details 
about this calculation in App. |X1 

In Fig. [3] we show the normalized density fluctuations for different modes as functions 
of X = kt], where 1] is the conformal time. The coupling between and 7'* results in 
damped oscillations as mentioned previously, and this is apparent in Fig. [3] for the modes 
that enter the horizon before kinetic decoupling. However, there is one major difference 



between this hidden sector scenario and the WIMP scenario studied in Refs. |52|, |5J] : the 
oscillation amplitude is much smaller than it is in the WIMP case. This traces back to 
the f^'j^ scattering cross section, which is ~ a|^/m^ ~ gw/''^wy ^^^^ 
weak scale, which is far smaller than the analogous visible sector cross section, which is the 
Thomson cross section ~ g^/f'^l- The smaller opacity in the hidden sector case results in 
large diffusion damping and hence the oscillations are highly suppressed. 

In Fig. m we plot transfer functions of the density perturbation of the hidden charged 
dark matter for various mx and combinations of (tan^(y, ^rh)- For smaller masses, kinetic 
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FIG. 3: The normalized amplitudes of dark matter fluctuation for different modes with comoving 
wavenumbers Xd = 0.01, 7.5, 75, 16 as functions of x = krj, where rj is the conformal time. We fix 
(tan6l^,^RH) = (V^,0.8) for this plot. 

decoupling happens later and the matter power spectrum is more suppressed. We define 
the cut-off wavenumber kcnt as the point where the transfer function first drops to 1/e of its 



value at small wavenumbers. For (tan^(y,^RH) = (y 3/5, 0.8) (top panel) and nix = 1 TeV, 
10 GeV and 100 MeV, we find fc^ut = 8.0 x 10^ Mpc'S 9100 Mpc-\ and 970 Mpc'S 
respectively. The free-streaming damping scale is 



where rjd is the comoving horizon at the time of kinetic decoupling, and rj^q is the comoving 
horizon at matter-radiation equality, and F(l/2) = y/n. For mx = 1 TeV, 10 GeV, and 
100 MeV, we find Ag^ = 9.3 x 10^ Mpc-\ 3.8 x 10^ Mpc~\ and 1500 Mpc-\ respectively. 
We see that for these masses, Afs < k~^^ — the cut-off scale is determined by the acoustic 
damping. 

In Fig. m we also show the transfer function of the dark matter density perturbation for 



the cases (tan^(jj/, ^rh) = (y 3/5, 0.1) (middle panel), and (10,0.1) (bottom panel). Com- 

paring with the top panel, we can see the transfer function cuts off at larger wavenumber for 
colder hidden sectors. Kinetic decoupling happens earlier and only the small scale modes 
that entered the horizon before kinetic decoupling get suppressed. For mx = 1 TeV, the 
cut-off is determined by the acoustic damping scale. However, for mx = 10 GeV and 100 
MeV, the free-streaming damping scale is comparable to the acoustic damping scale. 

The matter power spectrum contains all the information about linear Gaussian density 
fluctuations. These fluctuations are amplified by gravity to create non-linear structures, that 
is, dark matter halos. However, on scales below the cut-off the matter distribution is smooth 
and gravity cannot regenerate power on these scales. Thus the linear power spectrum also 
encodes information about the smallest building blocks of dark matter halos. The minimal 





(15) 




10 




lO'' 10" 
k[1/Mpc] 



-^10° 

E 

■8 



10" 



— 1 — I I I lllll 1 — I I I I nil 1 — I I I Mill 1 — I I I Mill 1 — I I I Mill 1 — I I I Mill 1 — I I I Mill 1 — I I I I III 



mx = 1 



(tane, ^rh) = (0.8, 0.1) 



m 1 — I I I Mill 

TeV 
lOGeV 
100 MeV 



A 



..I I I 



10° 10^ 10^ 10^ 10"* 10^ 10^ 10'' 10^ 



k[1/Mpc] 



lO'' 

E 

E 



^10-^ 



10 



-2 



—1 — I I I llll| 1 — I I I llll| 1 — I I I llll| 1 — I I I llll| 1 — I I I llll| 1 — I I I Mll| 1 I I I llll| 1 I I I llll| 

my = 1 TeV 



(tane, ^rh) = (10, 0.1) 



lOGeV 
100 MeV 




10° 10^ 10^ 10^ 10"^ 10^ 10^ 10^ 10^ 



k[1/Mpc] 



FIG. 4: Transfer functions of the normalized dark matter density perturbation amplitude for 
(tan6'(^,^RH) = (^3/5, 0.8) (top), (a/375,0.1) (middle), and (10,0.1) (bottom). 
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FIG. 5: Mass of the smallest virialized dark matter structure that can form, as a function of the 
dark matter mass mx for (tan(9|V,^RH) = (v/S/S, 0.8), (^3/5, 0.1), (10, 0.1), as indicated. 

mass of dark matter clumps may be estimated as 



where Vtm is the matter density today, and pcrit = SHq / (SnG) = S.lh"^ x 10"^'' GeV^ is the 
present critical density. 

In Fig. [5] we show the characteristic mass of the smallest objects as a function of dark 
matter mass mx for various values of {tan6'^,C,KR). For the lower reheating temperature 
parameter ^rh, ^cut is smaller by a few orders of magnitude for a given dark matter mass. 
For cold hidden sectors, kinetic decoupling occurs earlier, corresponding to smaller comov- 
ing horizons and larger fccut, which leads to smaller Mcut. We can understand this more 
quantitatively. The parameter kcnt is related to the comoving horizon when kinetic decou- 
pling happens as kcutVd = Xd, where a typical value is Xd ~ 5.9 in our case. On the other 
hand, we have rjd ~ l/^kd, so kcnt is proportional to the decoupling temperature, that is, 
kcvLt ~ ^kd, and therefore M^ut ~ (^kd)~^- As an example, consider the mx = 1 TeV case. 



From Fig. [21 the ratio of visible decoupling temperatures for the (tan6'(^, ^rh) = (v 3/5, 0.8) 



and (y3/5,0.1) cases is 7 x 10 ^, and so the expected ratio of minimum halo masses is 

3 X 10^ as may be ascertained from Fig. [51 

Our results in Fig. [5] show that the minimum halo mass — the building blocks of structure 
in the universe — could range anywhere from 10~^°Mq to the size of the smallest galaxies 
observed in the universe for dark matter masses within the MeV to TeV range. We will see 
in Sec. IVIII that if we restrict our attention to tan9^r < 10, we will require mx ^ 1 GeV to 
obtain the right relic density and be consistent with the observed morphology of galactic dark 
matter halos. For this range of parameter space, we see from Fig. [SI that the minimum mass 
halos are less than about lO^M©. At the present time, there is no way to test for halo masses 
as small as this. These predictions for structure formation are therefore indistinguishable 
from those of WIMP models, and current observations do not place constraints on the hidden 
charged dark matter parameter space. 




(16) 
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In the future, it may become possible to probe the upper end of this range with strong 
gravitational lensing of radio-loud quasars [6l|. In addition, if the dark sector is linked to 
the visible sector through connector particles with both visible and hidden sector quantum 



numbers [12|, ll3|, [ij], then self-annihilations of dark matter within these mini-halos might 



provide a detectable signature for the Fermi gamma-ray experiment [6^ 



63|, m, m m, m 



V. BOUND STATE FORMATION AND SOMMERFELD-ENHANCED ANNIHI- 
LATIONS IN THE EARLY UNIVERSE 

Charged dark matter annihilation may be enhanced in two ways. First, at low velocities, 
its annihilation cross section is enhanced by the Sommerfeld effect |28|]; in Feynman dia- 
grammatic language, this enhancement arises from diagrams with additional U(1)em gauge 
bosons exchanged in the t-channel that are higher order in gauge coupling, but kinematically 
enhanced at low velocities 70, 7l|, T2\- Second, charged dark matter may form bound states. 



which then annihilate. These are independent effects that must be included separately. Both 
act to enhance annihilation, and one might worry that they have negative implications. For 
example, the resulting enhanced annihilation rates may reduce the thermal relic density 
after chemical freeze out to negligible levels or otherwise change the relic density in a way 
that is excluded by astrophysical observations. 

In full generahty, these effects may be evaluated by solving the coupled Boltzmann equa- 
tions for the number densities of the bound state B = {f^^f^~') and the dark matter particle 
X = f^, 



dt 
dnx 
~df 



+ 3Hnx = -nxFann - "-^rrec + ?^B^ion , (17) 



where 



r,ec = r(X+X- ^ B^^) = Cnx {a,,,v) 
1 B = 1 (i^ ^ 7 7 ) = — - — 

F,,, = F(X+X- 7^^) =CnxS {(Tav) , (18) 

are the recombination (bound state formation) rate, the bound state decay width, the ion- 
ization rate, and the Sommerfeld-enhanced dark matter annihilation rate, respectively. In 
Eq. f|T8l) . C = (p^)/(p)^ accounts for the clumping of dark matter into halos, cta is the 
tree-level S'-wave annihilation cross section, and 

S = ^1^2ih_ (19) 



is the Sommerfeld enhancement factor [28|, |73| , where v is the center-of-mass velocity of the 



incoming particles. For v > ax-, 5" ~ 1, but for v <^ ax, S ~ nax/v may be a significant 
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enhancement. The recombination and photoionization and cross sections are related through 
the Milne relation 

where gf = 1 and gn = ^ are the statistical weights of the staus and the n-th level of the 
bound state, and i?„ = a^mx/4n^ is the binding energy of bound state level n. Photoion- 
ization may, in principle destroy bound states before they annihilate. 

The recombination rate requires some analysis. First, note that a photon is radiated 
in the recombination process; bound state formation followed by bound state decay yields 
a different final state than Sommerfeld-enhanced annihilation, and these two effects must 
be included separately; as we will see, both yield an enhancement that is parametrically 
~ ax I V for small velocities. 

The recombination cross section may be determined by using a hydrogen atom wave- 
function solution in the Schrodinger equation for a particle moving in a central potential 
V{r) ~ —ax jr. The recombination cross section for the n-th shell is [tJ] 



2«7r^ ax I B., ( V e__ 

3 mxmxmxv^\Bn + \mxv^ -2-,/^ ' ^ ' 

This result includes a correction factor from the fact that the electromagnetic force is long 
range, so that the incoming state is not a plane wave. The phase space distribution of the 
cold dark matter is given by solving the Boltzmann equation of Eq. ([6]). The solution is 
a Boltzmann distribution with dark matter temperature Tx ~ before kinetic decoupling, 
and an effective temperature 

T.^^-^ (22) 

after kinetic decoupling. The thus has Maxwellian velocity distribution at the effective 
temperature Tx, 



/ ™ \ 3/2 

f{v)dv = Ani^-^-^j v^e~^^dv. (23) 

Summing over all binding states, the thermally averaged total recombination cross section 
is 



rec ^ 



^ 4 2V 1 ax 
^ 3 vim\ 




(24) 



where we have made the replacement fg = ITxIvtix- We see that in the small kinetic 
energy limit, the thermally averaged recombination cross section is the thermally averaged 
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FIG. 6: The recombination rates (bound state formation rates) Tj-^c (solid) in the early universe as 
functions of the visible (CMB) temperature for dark matter mass mx = 1 GeV and (tan 0^^, ^rh) = 
(10,0.1) (upper) and mx = 1 TeV and (tan0(y,^RH) = (\/3/5, 0.8) (lower). Also shown are the 
Sommerfeld-enhanced annihilation rates Fann (dotted) and the Hubble expansion rate H (dashed). 



annihilation cross section of Eqs. (jlj) and ([5]), enhanced by a "bound state enhancement 
factor," that is, {a^ecv) = S^-cc {ctav), with 

Srec ~ 0.7 {-^) f ^) , vo/ax « 1 . (25) 
V87rao/ \ vq J 



The numerical prefactor is, for example, 0.5 and 0.05 for taxiO^ = y 3/5 and 10, respectively. 
Summing to higher binding states would give a factor of ~ 1.6 at most. 

In Fig. [6] we plot the bound state formation rates Free for mx = 1 GeV and 1 TeV. After 
chemical freeze out and before the bound state formation can take place, the cold dark 
matter number density is nx = ^x PcritO'~^ /f^x oc T^. Here we compare in the era before 
matter-radiation equality where the dark matter particles have not clumped into halos yet. 
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Next we will consider enhanced annihilation and bound state formation in the earliest bound 
objects, the protohalos. The bound state formation rate for different dark matter masses 
mx may be understood as follows. In general, there are three stages: 

1. Before kinetic decoupling the dark matter most probable speed is Vq oc T^/^. When 
Vo/ax ^ 1 the thermally averaged bound state formation cross section goes as 
(o'vecv) oc 1/fQ and hence the rate T^ec oc T^/^. 

2. After kinetic decoupling, the dark matter most probable speed becomes Vq oc T. When 
Vo/ax ^ 1 still holds, the bound state formation rate is Free oc Vq^ = constant. 

3. When the dark matter velocity drops to the point Vo/ax ~ 1, (crecv) becomes oc I/vq 
(cf. Eq. (121]) )• As a consequence Free decreases as T^, the same as the Hubble rate in 
the radiation-domination era. 

Also shown in Fig. E] are the Hubble expansion rate H and the Sommerfeld-enhanced 
annihilation rates Fann- The results can be summarized as 

H > Fann ^ F,ec (26) 

for mx = 1 GeV — 10 TeV. The bound state decay rate F^ and photoionization rate Fion 
can be much larger than the Hubble rate. But since bound states rarely form, we may set 
ub = in Eq. f|T71) . We see that for v -C ax, bound state formation has qualitatively 
the same effect on annihilation as the Sommerfeld effect — they enhance annihilation by 
factors of S and Srec, which are both proportional to ax/v, with order one coefficients. Since 
the Sommerfeld enhancement factor is also valid for v ^ ax, while the bound state effect 
becomes negligible at large v, one may draw qualitatively correct conclusions by considering 
only the Sommerfeld effect and neglecting the effects of bound state formation, but we 
emphasize that these effects are independent, yielding two different final states. 

We may now answer the questions posed at the beginning of this section. Because the 
Hubble expansion rate is larger than the Sommerfeld-enhanced and bound state-catalyzed 
annihilation rates, the thermal relic density is not modified significantly between the times 
of chemical freeze out and matter-radiation equality. 



VI. SOMMERFELD-ENHANCED ANNIHILATIONS IN PROTOHALOS 

We now consider these non-perturbative effects on the dark matter annihilation after 
matter-radiation equality, when the growth of structure enters the non-linear regime. The 
smallest structures undergo gravitational collapse first and the more massive ones form later. 
The minimum halo mass Mcut is set by the cold dark matter kinetic decoupling temperature 
and we have computed this for the hidden charged dark matter in Sec. |IV] (cf. Fig. |5]). 
The redshift at which these objects form (virialize) depends on the underlying matter power 
spectrum jTsj, which we have also computed in Sec. [IVl Simulations 7g| find Zc ~ 30 for 



a typical 10~^ Mq halo and we will assume this value in this work, but note that in detail 
the collapse redshift depends on Mcut- We take the overdensity of the virialized region to 
be ~ 178 times the ambient cosmological density at (appropriate for WMAP cosmology) 
and estimate the velocity dispersion of the dark matter particles in the halo as: 

1/3 



.o = \/^-3.0xlO f^j (1 + ..)/ . (27) 
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Hidden charged dark matter particles in our model with masses mx ^ 1 GeV have velocity 
dispersion vq -C ax in these protohalos. They thus form bound states and pair-annihilate 
with an enhanced cross section, from an initial number density nx — 6 (mx/GeV)~^ cm~^ 
for about a Hubble time before the protohalos merge into larger halos. In these larger halos, 
the dark matter particles have smaller phase space density and so the annihilation rate is 
lower. 



Dark matter annihilation in protohalos was used in Ref. [32| to obtain stringent con- 
straints on models with the Sommerfeld enhancement effect. In that work dark matter 
particles were assumed to annihilate dominantly to visible photons and e^e~ pairs, which 
contribute to the extragalactic diffuse gamma-ray flux and generate CMB temperature and 
polarization power spectra. However, in the present context, where the hidden sector is not 
coupled to the visible sector, only a much weaker, separate, constraint from CMB anisotropy 
measurements applies. The CMB measurement of the matter density VL^h'^ can be compared 
with those obtained by distance measurements using Supernovae Type la (SNIa) and baryon 
acoustic oscillations (see, for example, Ref. 77|). A rough estimate will suffice for our pur- 
poses here and we will assume that no more than 10% of the total matter can be converted 
into radiation between the epoch of recombination and now. 

This constraint may then be applied to our models, where the S'-wave stau annihilation 
cross section of Eq. (jlj) is enhanced by the Sommerfeld factor S. As mentioned before, when 
vq/o.x -C 1 the bound state formation cross section is of comparable size. One demands the 
fraction of stau annihilation in the protohalos 

f ^nxS{aAv)Ct^ 178^^cPcrit + ^^)% (^^ ^) < 0.1, (28) 

mx 

during t ~ 5.6 x 10^^ (1 + z^~^-^ s, the age of the universe at redshift Zc ^ 1. The density 



profile of the minimum mass halos [76| is expected to be similar to those of the larger 
halos j75|, and we include an appropriate clumping factor C = 7 for the minimum halo mass. 
We find that for {tan 9^, ^rh) = (10, 0.1), / ~ 10"^ - 10"^ for mx = 1 GeV - 10 TeV and 



ax fixed by Qxh"^ = 0.11. For (tan^jV.^Rn) = (^3/5, 0.8) and (^3/5,0.1), / ~ 10"^ -10"^ 
for the same mass range and the corresponding ax's. These values are far below the bound 
of / < 0.1, and so the annihilation of dark matter in protohalos does not constrain the 
scenario in the case of purely hidden charged dark matter. However, if there were connector 
particles with both hidden and visible sector quantum numbers mediatiiig annihilation to 
observable particles, a much more stringent bound / ^ 10~^ may apply [32] • This is from 
the consideration that (visible) photons ejected at an energy > 300 GeV or < 100 keV 
should not lead to heating and ionization of the intergalactic medium that contradicts CMB 
and large scale structure observational data. 

More massive halos have larger dark matter velocity dispersions. For example. Milky 
Way-size halos have vo ~ 270 km/s ~ 10~^, while dwarf galaxies can have velocities as 
low as vq ~ 10 km/s ~ 7 x 10~^. For hidden charged dark matter particles with masses 
mx = 1 GeV — 10 TeV, then, only a fraction / ~ 10^^ annihilates during the age of the 
Universe. 



VII. ELLIPTICITY AND CORES OF DARK MATTER HALOS 

At late times, elastic scattering between dark matter particles in the halos of galaxies 
or clusters of galaxies may change the shape of const ant- density contours and lead to the 



17 



formation of a core, that is, a central region with constant density. The main effect results 
from transfer of kinetic energy in collisions and if this process is fast enough to create 0{1) 
changes to the energy of the dark matter particles in the halo, then it will drive the halo 
towards isothermality and isotropize the mass distribution. 

The true dynamical picture is more complicated. Initially, heat is conducted from the 
hotter outer to the cooler inner parts of the halo through collisions. This heats up the core 
causing it to puff up, which gives rise to a density profile that is much shallower (or even fiat) 
compared to the initial cuspy density profile of the central regions of the halo. In addition, 
self-annihilations will also lead to cores [tI], however our calculation in Sec. |Vl]has shown 
that this is a not large effect. The collisions also erase velocity correlations and lead to more 
a more spherical (rounder) core. We find that the observed ellipticity of galactic dark matter 
halos provides the strongest constraints on charged dark matter models. 

Over periods long compared to the relaxation time, ejection of dark matter particles 
from the parent halo in collisions will further cool the core. For an isolated halo, this will 
eventually lead to core collapse and result in a steep cusp. In a cosmological setting, however, 
halos accrete dark matter particles and this would offset the energy loss due to collisions. 
Cosmological simulations [tq], S^I seem to indicate that core collapse occurs on a time scale 
much longer than the relaxation time scale, if at all. We will, therefore, not consider this 
effect any further here. 

Shapes of dark matter halos of elliptical galaxies 8l|, l82| and clusters SJ, l83| are decidedly 
elliptical and this fact may be used to put constraints on the self-interactions. Such an 
analysis [3^] was carried out in the context of the self-interacting dark matter proposal js^] . 
This proposal was designed to explain why observations of the rotation curve of low surface 
brightness galaxies seemed to indicate that the dark matter distribution had a fiat density 

' I, 83] and to explain the observed census of dwarf galaxies in 
The status of these discrepancies is unclear — while it seems clear 



profile in the inner core 
the Local Group [ 



85 



89 



that the census of the observed local group galaxies are broadly in agreement with the cold 
dark matter model [OO, W], the fit to rotation curves of low-surface brightness galaxies is 



still problematic |92|, ,93,, \M M 



The strongly self-interacting dark matter proposal motivated the first nu merical sirnula- 

loi, 



33, 79,^ 



99, 100, 101 



tions to deduce the effects of dark matter self-interaction ^ 

T03l and they validated the expectation that during the evolution to isothermality, the dark 
matter core becomes rounder. These simulations indicated that only regions of parameters 
space with 500 GeV~^ < o-Du/mx ~ 5000 GeV~^ (0.1 cm^/g < (JDM/mx ^ 1 cm^/g) could 
introduce observable features on the scale of the dwarf galaxies while at the same time be- 
ing consistent with observations of larger galaxies. In the constraint quoted above, ctdm is 
the dark matter elastic self-scattering cross-section modeled as a hard-sphere interaction. 
The analysis of the shape of the dark matter halo of a particular cluster of galaxies (using 
gravitational lensing) indicated that much of the above preferred region was ruled out j33 |. 

We revisit the constraint from the inferred ellipticity of dark matter halos in the rest of this 
section. These constrain t s ar ise from a wide-ran ge of observati ons including X-rays 8l|, l83| , 

Recent wor k on combining 
in clusters 



strong lensing 104, 105, 106l | and weak lensing JiOTI . 108, 109 



110 



as 



different measurements to reveal the anisotropy of velocity dispersion 
opposed to the shape of dark matter distribution) is a different but complementary way to 
constrain the self-interaction of dark matter. 

To estimate how these observations may be used to constrain the dark sector Coulomb 
interactions, we calculate the relaxation time for establishing an isothermal halo. We will 
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then assume that the time scale for isotropizing the spatial distribution of the dark matter 
halo is the same as this relaxation time and use constraints from measurements of the 
ellipticity of galaxy halos to put limits on the dark sector Coulomb interaction. The notable 
feature of the Coulomb interaction is the strong dependence of the energy transfer rate on 
the relative velocity of the interacting particles. This ties in with other investiga tions that 
considered ve l ocity -dependent interaction cross sections. See, for example, Refs. 111 . 112l . 
TT3l.ITli.fTTi.ITTi. 



The Rutherford scattering cross section is 



da a\ /„ V 

^ (29) 



where 6 is the scattering angle in the lab frame. We assume the dark matter particles in the 
halos have a local Maxwellian velocity distribution f{v) (normalized to unity) and velocity 
dispersion (u^) = (3/2)fQ(r), which in general varies spatially within a halo. The kinetic 
energy exchange in each colhsion is SEk = Ek {I — cos 9), where Ek = mxv'^/2. The rate of 
energy transfer is, then, 

Ek = j dvdil^f{v)6EkVnx 

a\r 4 11 Px f , , /il — cos 9 „2 /„2 
= 27T-^^^-mx— dvdcos9 — , , , t;e~" 
Am^ y/^ vi 2 mx J sin^ (9/2) 

- ln(l - cos6'min)t^o • (30) 



mx Vq 

The minimum scattering angle is related to the maximum impact parameter through 

CXx 

&max = 2 (^min/2) , (Sl) 

nixVQ 

where 6max should be chosen to be 

mx Vq . . 

V47r ax px 

the Debye screening length in the f'^ plasma. We find constraints on ax by demanding that 
the relaxation time be larger than the age of the universe, 

r, = Ej,/Eu ~ , r-\ 1^ ^ 1 > 10^° years , (33) 




A^a\ px 

where the "Coulomb logarithm" is ~ 90. The constraints on ax obtained above scale 
inversely with the square root of the "phase space density" Q = Px/vq, and this indicates 
that the best constraints may be obtained by studying the cores of galaxies rather than 
clusters of galaxies. 



Many elliptical galaxies show clear evidence for flattened, triaxial dark matter halos [81 
Using the profi l es of the total mass enclosed M(r) and the halo concentration parameters c 
from Refs. |ll7 . 118l |. we derive the radial velocity dispersion f ^(r) = (3/2)i;g(r) and the dark 



matter density px{f) at a radius ~ 3 — 10 kpc. The dark matter density drops from 3.5 to 0.7 
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GeV/cm^ as one moves outwards, while the velocity dispersion decreases slightly from 270 
to 250 km/s. Using these values in Eq. (!33l) . we obtain a very stringent bound on ax- This 
constraint is shown as the lower solid line in Fig. [1] If we demand that the WIMPless scenario 



provide the right rehc abundance of — 0.11, then for (tan^{y,^RH) = (y3/5,0.8) and 

(y^3/5, 0.1), nix ~ 100 GeV. However, for (tan6'^,^RH) = (10,0.1), hidden charged dark 
matter as light as rrix ~ 1 GeV may freeze out with the correct thermal relic density without 
being in conflict with the constraint from elliptical halos. 

The constraints obtained above are comparable to those deduced in Ref. ^B]- This 
previous work demanded that the properties of the Milky Way's dark matter halo (and by 
extension that of other galaxies as massive as the Milky Way) should not deviate by order 
unity from those predicted by coUisionless dark matter simulations. If we compare their 
result to our Eq. fl33|) with px = 0.3 GeV cm~^ and vq = 200 km s~^, we obtain weaker 
constraints on ax by a factor of ~ 6 at the same mx- We trace the bulk of this discrepancy 
to the fact that they assume total dark matter mass of IO^^Mq in their calculation and 
neglect the contribution of stars in setting the local dark matter velocity dispersion. The 
constraints we have obtained from detailed observations of an elliptical galaxy with strong 
X-ray emission resolved down to 3 kpc of the center — better than we can do within our own 
galaxy for now — are comparable to the bounds quoted in Ref. [l9|. Future o bserv ations of 
hyper- velocity stars within the Milky Way could strengthen these constraints liol . 



We may also consider dark matter halos more massive and less massive than that of the 
Milky Way. Cluster halo shapes are measured with X-rays and strong gravitational lensing 
at radii ~ 100 kpc. We consider some Abell clusters for which the ellipticity profiles are 



determined in Ref. [83|. The radial velocity dispersion at this scale is typically \ Kv^) 



1000 km/s, as inferred from the M(r) profiles determined in Ref. [l20l |. As a consequence 
the bounds on ax are about two orders of magnitude weaker. 

In principle, the most stringent constraint using Eq. (!33|) would be from smaller spiral 
galaxies and the Local Group dwarf galaxies, where we expect larger dark matter phase space 
densities and therefore shorter relaxati on t imes. Current observations of the shapes of dark 
matter halos of nearby spiral galaxies jl21 | find, however, that their gravitational potential 
is quite round, and that they seem to prefer a central core. It is interesting to note that 
within our model, the effect of dark sector Coulomb interactions could leave a dynamical 
imprint in some of these nearby spiral galaxies but not the larger elliptical galaxies. We 
leave more detailed investigations for future work. 

The case of the dwarf galaxies in the Local Group, that is, the satellites of the Milky 
Way and Andromeda galaxies, is particularly interesting. The stars in these galaxies have 
very small velocity dispersions — of the order of 10 km/s or even smaller. Combining 
this with the observed extent of the stellar population, one may infer the mass of the dark 
matter within the stellar extent of these dwarf galaxies. The results of such an analysis show 
that the dwarf satellites of the Milky Way are consistent with dark matter central densities 
of about O.IMq/pc^ ^ 4 GeV/cm^ jl22 |. Unlike the galaxies we have been considering 



previously, these dwarfs present a complication. Their present day properties are set rather 
dramatically by the Milky Way galaxy. As they fall into the Milky Way, gravitational tidal 
forces will strip them of mass on the outside and thus reduce the velocity dispersion of the 
dark matter particles inside the satellite's halo. Exactly what the dispersion is depends on 
the initial mass, the extent of the tidal mass loss and dark matter self-interactions. If the 
dispersion is of the order of 10 km/s (similar to the observed stellar velocity dispersion), 
the energy transfer time scale is short in the dwarf satellites for ax and nix that are at the 
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edge of the allowed region in Fig. [H For these values of mx and ax , interactions will almost 
certainly form cores in the dwarf galaxies . S u ch c ores are compatible with observations 



unless they are larger than about 300 pc [l23l . Il24l |. Future astrometric measurements of 



individual stars in these satellite galaxies will be able to measure the density profile of dark 



matter halos on these scales 125 



Surprisingly, if we use the relation between the core radius and core density for a self- 
gravitating i sothe rmal sphere, pcore = 9'f^o/('^^^'^core) with vq = 10 km/s and pcore — 
4 GeV/cm^ jl22l |. we obtain a core size of 200 pc. However, our estimate has not ac- 
counted for some crucial factors. If the mean free path is much shorter than the typical 
orbit of the dark matter particles within the satellite, heat conduction will be suppressed. 
In addition, the passage through the Milky Way will introduce interactions with the dark 
matter in the Milky Way and this could change the density profile of the satellite. 

More detailed work is required to consider the effect of dark sector Coulomb interactions 
on the satellites of the Milky Way and other galaxies. We conjecture here that there would 
be regions of allowed parameter space where the satellites galaxies and other small field 
galaxies would show constant density cores in the center and reduced substructure within 
their halos. This is beyond the scope of the present work, but we urge the reader to keep in 
mind that future work in this direction could lead to interesting astrophysical phenomena 
and perhaps also rule out some of the allowed parameter space. 



VIII. DARK MATTER IN THE BULLET CLUSTER 

The Bullet Cluster is a rare system where a subcluster is seen to be moving through a 
larger cluster. Through observations in the optical and X-ray and strong and weak gravita- 
tional lensing observations, astronomers have been able to map out the spatial distributions 
of the stars, gas and dark matter in this system. From these inferred distributions, it is 
clear that dark matter tracks the behavior of stars, which are coUisionless, rather than 
the gas. This observation has allowed stringent bounds to be placed on the dark matter 
self-interaction strength. 

With respect to the self-interacting dark matter proposal js^] discussed in the previous 
section (velocity- independent cross section), the Bullet Cluster observations demand [s^, 
1261 that aBu/'mx ~ 3000 GeV ^ {o'DM/'mx ~ 0.7 cm^/g). These are the most direct 



constraints on the self-interaction of dark matter. 

These bounds have been derived based on different considerations, including the observed 
gas and dark matter offset, the high measured subcluster ve locit y, and the survival of the 
subcluster after having moved through the Bullet Cluster |l26 |. The last phenomenon. 



survival of the subcluster, turns out to provide the s trongest constraint and hence we will 
focus on that. We follow the approach of Ref. |l26 | to derive bounds on the dark sector 



Coulomb interaction, but relax the assumption of a hard-sphere interaction cross section. 

The subcluster experiences a net particle loss in a collision when the velocities of both 
particles in the main and subcluster become larger than the escape velocity f esc — 1200 km/ s. 
This condition may be turned into an effective range for the scattering angle in the subclus- 
ter's reference frame 

^<cos0<fl-^^ , (34) 
vi \ n J 

where Vi ~ 4800 km/s is the velocity of the main cluster incoming particle before the 
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collision. We now assume that the subcluster passed through the main cluster's center so 
that it saw a surface density ~ 0.3 g/cm^. Demanding that the fraction of particle loss 
be no greater than 30%, we have an upper bound on ax'- 




dn — 



da 2tt a\ 



1 



1 



(35) 



dQ mx Tn\ vf 



1 - cos 6'min 



1 — COS 6, 



< 0.3 . 



max J 



This bound is given in Fig. [H We see that it is about four orders of magnitude weaker than 
that obtained from considerations of the ellipticity of galactic dark matter halos. Equation 
( !35|) suggests that an improvement of this kind of bound can be made by considering a galaxy 
falling into a cluster, which results in a larger surface density = / pxdl ~ 30 g/cm^, 
where the galaxy path is / ~ 27ri?tage/6rdyn- By demanding similarly that the galaxy does 
not lose more than 30% of its particles during orbiting in the cluster for about the age of 
the universe, we obtain a bound on ax roughly ten times stronger than from the Bullet 
Cluster, but still weaker than the bounds obtained from the observed ellipticity of galactic 
dark matter halos in the previous section. 

IX. CONCLUSIONS 

We have investigated the astrophysical and cosmological consequences of dark matter 
that is charged under an unbroken hidden U(1)em- We find that this is a viable and natural 
dark matter candidate. The salient arguments leading us to this conclusion are as follows. 

• We investigated the dark matter-dark photon and dark matter-dark neutrino scattering 
processes in the early universe and found that the dark matter kinetically decouples 
early and therefore behaves as cold dark matter. The predictions for the temperature 
at which dark matter kinetically freezes out are shown in Fig. O 

• Our calculations for the power spectrum of density fluctuations of dark matter showed 
that in these models (for mx > GeV) we expect structure all the way down to at 
least 10'^ Mq. The predictions for these building blocks of structure in the Universe 
are summarized in Fig. O At the present time, we have no way to distinguish these 
hidden sector charged dark matter models from canonical WIMP models using matter 
power spectrum observations. 

• Although the dark matter particles in our model are cold, they also have long-range 
Coulomb interactions (but with a smaller fine structure constant for masses below 1 
TeV). Since the hidden sector U(1)em is unbroken, we have equal numbers of positively 
and negatively charged particles and hence we investigated the formation of bound 
states. We found that for masses below 10 TeV, bound state formation and the 
subsequent annihilation does not change the relic density of dark matter particles 
appreciably. 

• The self-interactions mediated by the long-range (hidden sector) Coulomb force can 
affect non-linear structure formation, especially at the small mass end. We found 
that the most stringent constraints arise from the observed ellipticity of galactic dark 
matter halos. (The self-interactions, if too strong, will make the core of dark matter 
halos round.) 
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• Putting together all of the results above, we find that the mass of a dark matter 
particle charged under a hidden sector unbroken U(1)em must be larger than about 
a GeV, z/we restrict our attention to hidden sector weak mixing angles tan 6*^ < 10. 
Larger values of the weak mixing angle will allow smaller dark matter particle masses. 
The detailed constraints imposed by the galactic dark matter halo observations, the 
Bullet Cluster observations and the requirement of obtaining the right relic density 
are shown in Fig. [T] 

• We also found that smaller galaxies such as the observed low surface brightness spirals 
and satellites of the Milky Way and Andromeda galaxy are likely to have constant 
density cores in their dark matter distribution. More work is required to pin down the 
core sizes, but this is clearly the regime where large differences from the predictions 
of coUisionless cold dark matter, if any, will be found. 



One of the promising avenues forward that we are currently exploring is including a 
connector sector linking the hidden and visible sec tors. Such a connector sector will give 



rise to new phenomenology, further constraints 32l. ll27l | , and interesting implications for a 
variety of dark matter detection possibilities IsT 14. 17, 39]. 



A second promising avenue concerns structure formation on small scales. Our formalism 
suggests that in the allowed parameter space, the largest deviations from coUisionless cold 
dark matter will appear in galaxies much less massive than the Milky Way. A detailed 
investigation of this aspect of hidden sector dark matter models is beyond the scope of the 
present work, but simple estimates suggest that, for large regions of model space, small 
galaxies will form cores and the substructure in all dark matter halos will be reduced. These 
expectations are in stark contrast to the predictions of coUisionless cold dark matter models. 

We also stress that there are many other possibilities in the hidden charged dark mat- 
ter framework. We have considered the hidden stau as the dark matter. Other charged 
particles, for example, the hidden tau, chargino, quarks, and squarks, are also possible. 
This scenario also supports multi-component dark matter, with several hidden sectors, each 
with its own dark matter particle contributing significantly to the relic density through the 
WIMPless miracle. Alternatively, even with only one hidden sector, stable hidden staus and 
hidden (light, but massive) neutrinos could form mixed dark matter, with both hot and cold 
components. All of these scenarios merit further study. 

To summarize, we have investigated the astrophysics and cosmology of dark matter 
charged under a hidden unbroken U(1)em- We find that for dark matter masses larger 
than about a GeV, these models have the right relic density, make cosmological predictions 
currently indistinguishable from the usual WIMP models, and are consistent with obser- 
vations on galactic and cluster scales. At the same time, these dark matter candidates 
are coUisional and their annihilation is enhanced, implying consequences for future observa- 
tions and experiments that may differ radically from WIMPs and other more conventional 
possibilities. 
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APPENDIX A: CALCULATION OF DARK MATTER POWER SPECTRUM 

In this appendix, we outhne the calculation to compute the small-scale power spectrum of 
dark matter. In our hidden charged dark matter scenario, the relevant interactions are elastic 
scattering off both hidden sector photons and neutrinos. For simplicity of presentation, we 
present results below for just Compton scattering off of photons. This is the dominant 
process in most of parameter space. In our numerical analysis and the results presented 
above, however, we include also the neutrino scattering process. 

The cross section for Compton scattering is 



U /I Li jj^ / A 1 \ 

= ^ ■ ^^^^ 



This differs markedly from the usual WIMP interaction cross section with standard model 
fermions, which is suppressed relative to this by (T/m^)^. This would result in lower decou- 
pling temperatures although the effect is mitigated to a large extent by the smaller coupling 
required to obtain the correct relic abundance, as discussed in Sec. IIIII 

W e stu dy the Fourier-decomposed perturbations in the conformal Newtonian gauge. See 



Ref. |128| for details. For this purpose it is better to use the visible photon temperature 
T as the "clock." Different hidden sector reheating temperatures T^-^ = ^rh^rh result in 
very different dark matter kinetic decoupling temperatures T^d (cf. Fig. [2]). However, their 
effect on the relation of the conformal time r] = J dt/a{t) to the visible photon temperature 
is small. 

For simplicity we assume that the fluctuation of the gravitational potential is determined 
solely by the visible sector photons because the energy density of the visible sector dominates. 
With this approximation we may use the analytic result for the Newtonian potential in the 
radiation dominated regime. 



$ = 3$; 



sin(fc?7/v^) — (krj/y/S) cos{kri/\/3) 



(A2) 



where k is the Fourier mode of interest and $p is the primordial amplitude. This approxi- 
mation results in < 10% error, depending on ,^rh- The •y'^ perturbations may be expanded 
in Legendre polynomials and this decomposition results in a multipole hierarchy: 

Q'^ + kQ'l = -<i> 



^ \3 ^ 3 V 3 V 3^ 
^ 5 ^ 10 ^ 



©z = ^ [^0ti-(^ + l)0f+iJ -ref for/>3, (A3) 

where the 0f are the multipole moments of the hidden sector temperature field, to a good 
approximation = — $, f = anxcffh^h is the scattering rate, and nx is the number 
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density. The dark matter density and velocity perturbations are 



Sx + Ox 

Ox + - Ox 
a 



-3$ 



(A4) 



where R = 3px/4:p^h is the ratio of dark matter to hidden sector photon energy density. 



The interaction term 



3f 

'r 



45 



m 



(A5) 



X 



is the rate of transferring momentum \p\ ~ T'^ from the hidden sector photons to the dark 



matter particles (cf. the second term in Eq. (fT3|) ). We follow the treatment in Ref. S^] to 
set the dark matter sound speed Cs and shear ax- 

The main difference between our hidden sector scenario and standard cosmology is that 
the hidden sector photon decoupling epoch is much earlier than the drag epoch of the hidden 
charged dark matter. We include hidden sector photon higher multipoles and truncate the 
series at /max = 10 to accurately account for the free-streaming of photons. 

To obtain the transfer function we solve the dark matter fluid equation well into the 
decoupling regime and the evolve further in time (when the calculation becomes computa- 
tionally intensive) using the free-streaming solution j52| 



Sxiv) = exp 




A I ^ 




(A6) 



where kj = rjlT [r]^) /tux, to the time of matter-radiation equality. 
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